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At extreme pressures and temperatures, such as those inside planets and stars, common 
materials form new dense phases with compacted atomic arrangements and unusual physical 
properties. The synthesis and study of new phases of matter at pressures above 100 GPa and 
temperatures above 10 4 K— warm dense matter— may reveal the functional details of planet 
and star interiors, and may lead to materials with extraordinary properties. Many phases 
have been predicted theoretically that may be realized once appropriate formation conditions 
are found. Here we report the synthesis of a superdense stable phase of body-centred- 
cubic aluminium, predicted by first-principles theories to exist at pressures above 380 GPa. 
The superdense Al phase was synthesized in the non-equilibrium conditions of an ultrafast 
laser-induced microexplosion confined inside sapphire (a-AI 2 0 3 ). Confined microexplosions 
offer a strategy to create and recover high-density polymorphs, and a simple method for 
tabletop study of warm dense matter. 
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Ordinary substances behave in extraordinary ways at high 
temperature and pressure, resulting in completely new 
phases of matter 16 . This type of matter regulates processes 
deep in the Earth as well as in the cores of other planets 1,7,8 . The quest 
to simulate these conditions, and study the behaviour of materials 
under extreme pressure and temperature has a long history 9 , with 
a dramatic advance in this direction being made by the invention 
of diamond anvil cell 10 . During the last century, pressure and tem- 
perature records have increased continuously, and it is now possible 
to obtain multimegabar static pressures, with simultaneous heating 
by laser up to 3,000 K (refs 1,5,11,12). In fact, recently, conditions 
at the Earths core of 364 GPa and 5,500 °C were recreated in the 
lab 13 . Even more extreme but transient high-pressure/-temperature 
conditions have been created by dynamic methods using explosives 
and megajoule- class lasers 1419 . 

It has been predicted that new classes of materials with 
unusual combinations of physical properties should exist under 
conditions of extreme pressure and temperature 4,20 22 . For instance, 
high-pressure structural behaviour of crystalline metallic Al phases 
has been the subject of first-principles theoretical studies 20,23,24 , 
which predict a sequence of phase transitions with estimated static 
transition pressures within the 120-360 GPa and 200-560 GPa 
ranges for fcc=>hcp (face-centred-cubic to hexagonal-close-packed) 
and hcp=>bcc transitions, respectively. However, only the exist- 
ence of hcp-Al phase has been confirmed experimentally so far 22,25 . 
Prediction of aluminium phases at terapascal pressures has been 
reported 26 but could not be validated by state-of-the-art high-pres- 
sure synthesis methods. 

Recently, it was demonstrated that it is possible to create extreme 
pressure and temperature conditions in table-top laboratory experi- 
ments with ultrashort laser pulses focussed inside a transparent 
material, with intensity in the focal spot significantly above the 
threshold for optical breakdown 27 29 . Tight focusing of a ~100-fs 
laser pulse with 100 nj energy deep inside a crystal, deposits an 
energy density of several MJcm -3 into a submicron volume, several 
times higher than the strength of any material. This causes the solid 
to superheat and form a plasma. The confined plasma explodes 
and generates a powerful shock wave that expands out of the focal 
volume and compresses the surrounding pristine material 27 29 . The 
shock wave of the expanding plasma can reach pressures of over 
10 megabar (1 TPa), triggering dramatic changes in the material 
properties. One of the essential novelties of the new method is that 
the modified material remains confined in a strongly localized region 
inside a bulk of pristine crystal, and can be investigated later. 

Here we report on the synthesis of a new stable high-pressure 
phase, bcc-Al, formed by ultrafast laser-induced microexplosion 
confined inside sapphire a-Al 2 0 3 (refs 28,29). The phase is pre- 
served inside the amorphous compressed sapphire produced from 
a hot, dense, non- equilibrium and short-lived plasma. We suggest 
that the sapphire constituents, aluminium and oxygen, become 
spatially separated in the hot expanding plasma because of the 
differences in the ion diffusion velocity, and a new bcc-Al phase 
is 'frozen within the compressed sapphire. Synchrotron X-ray 
diffraction microanalysis (/iXRD) shows that the diffraction peaks 
observed in our sample perfectly match the bcc-Al structures 
predicted to exist at high pressures 20,23,24 . 

Result 

Formation of voids. In our experiments, we irradiated a-Al 2 0 3 
(corundum, a variety of sapphire) crystals with 150-fs laser pulses 
focused at depths of several micrometers, using a high numerical 
aperture microscope lens (Methods). Sapphire has a very large 
bandgap of 9.8 eV, and is of technological importance, given that 
it is a high-hardness ceramic material with a high bulk modulus 
of B = 239-257 GPa (ref. 30), and because the void inside 
sapphire is well distinguishable by an optical contrast and well 




Figure 1 1 The method to create high-density phase by the ultrashort 
laser pulse triggered and spatially confined microexplosion. A tightly 
focused femtosecond pulse generates hot plasma in sapphire, which in 
turn produces a strong shock wave, which compresses the material against 
the surrounding pristine Al 2 0 3 crystal into a densified amorphous phase A 
and forms a void V in the centre of a laser-induced microexplosion. 
The superdense bcc-Al phase was found within this compressed material. 



defined when compared with softer materials such as silica or 
polystyrene. 

Laser radiation caused the formation of voids with compressed 
shells within the sapphire crystal (Fig. 1). The size of the voids and 
the compressed- shell thickness were measured using transmission 
and scanning electron microscopes (TEM & SEM) after opening 
the enclosed laser- affected zone by Focused Ion-beam milling. 
The average radius of the combined void and shell structures 
was -360 nm, twice as large as the shell thickness of ~180nm. 
Assuming conservation of mass, the average shell density was 
5=1.14 p s = 4.54gcm~ 3 , where p s = 3.98gcm~ 3 is the initial density 
of sapphire. This shell was identified as amorphous material by TEM 
and chemical etching because amorphous sapphire has a much 
higher solubility in hydrofluoric acid than crystalline sapphire 28 . 

Samples for /iXRD analysis were prepared by making regular 3D 
arrays of laser-modified areas with a period of 2 jam within a region 
of (50x50x70) Jim 3 in polished slabs of a-A\ 2 0 3 . ^uXRD characteri- 
zation of the laser-processed region was carried out using a mono- 
chromatic 33keV (A = 0.3757 A) X-ray beam with a cross-section of 
5x10 [im at full width half maximum (FWHM). The microfocused 
X-ray beam was oriented axially, approximately parallel to both 
the original laser beam and to the crystallographic z-axis of the 
sapphire, and was scanned in 10 Jim steps over both patterned and 
non-patterned regions. 

Evidence of a superdense phase of Al. Figure 2a,b show 2D ||XRD 
images collected at laser-modified and intact areas in the sample. 
There are pronounced diffraction rings in the modified regions, 
whereas pristine regions show only a background of single-crys- 
tal-like sharp spots. An azimuthally-integrated |lXRD pattern of 
Figure 2a was thoroughly compared with the X-ray Bragg peak 
positions of known and theoretically predicted high-pressure A1 2 0 3 
phases 31 , and for all possible Al-O compounds at various combi- 
nations of density and temperature, but no satisfactory match was 
obtained. Also, given the large difference in electronegativity at high 
pressure 32 and in the atomic radii 33 of Al and O, a solid solution of 
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Figure 2 | Evidence of bcc-AI formation within the laser structured 
regions. Micro-X-ray diffraction (jiXRD) image acquired (a) in the centre 
of the Shockwave compressed area, and (b) outside the Shockwave 
compressed area, (c) Comparison between a radially integrated 
experimental fiXRD profile obtained from the experimental data shown 
in (a), and theoretical ,uXRD profile expected for bcc-AI refined using 
Materials Studio package, assuming 18 nm crystalline size. The crystalline 
(hkl) plane indices are shown next to the peaks. For comparison, simulated 
profiles of host sapphire (vertically offseted grey) and native fcc-AI 
(purple) are also shown. 
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Figure 3 | Rietveld refinement profile of the bcc-AI structure. A difference 
curve is plotted for comparison with offset. The last Bragg peak intensity 
was not included in the fit, because only a fraction of the diffraction ring 
was acquired in the XRD image (enclosed in a shaded-box). 



Al and O cannot explain the observed /iXRD spectra, as may be 
seen from the empirical Hume-Rothery analysis 34 outlined in 
the Methods. Further examination revealed that the measured 
peak positions could be perfectly indexed as a bcc-AI structure 
(space group Im-3m) with a lattice constant a = 2.866 A (11.77 A 3 
per atom), which matches a predicted 20,23,24 , but not yet experi- 
mentally observed, high-pressure Al phase; this bcc-phase is 41% 
more dense than the native fcc-phase of Al; the density of bcc-AI 
p bcc = 1.41 xp fcc = 3.81 gem -3 . Figure 2c shows a comparison between 
the experimental /iXRD pattern and Rietveld- refined profile of 
bcc-AI. The refinement fit is shown together with experimental in 
Figure 3. Shapes of the Bragg peaks were fitted using Rietveld refine- 
ment with crystallite size and microstrain as fitting parameters, 
yielding a crystallite size of 17.9 nm and strain of 0.26%, closely 
matching the values calculated using Scherrers algorithm 35 : 
L = KX/(j5 cos 0), where L= 18±2nm is the size of the nanocrystal- 
lite, K= 0.891 is a constant, /3 is the FWHM of the diffraction peak, 
and 6 is the diffraction angle. 
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Figure 4 | Determination of the residual pressure in the vicinity of the 
retrieved bec- and fcc-AI nanocrystallites. (a) Micro-X-ray diffraction 
(jiiXRD) patterns from the laser-structured region (the pattern of Figure 2) 
where bcc-AI signatures or fcc-AI are dominant, (b) Comparison of the two 
integrated ^uXRD patterns shown in (a). Numerical fit of the experimentally 
observed peaks by the fcc-AI structure with a cell size of a = 3.584 A is 
achieved; bottom purple marks show the expected position of the fcc-AI 
peaks at room conditions with a = 4.049 A. The retrieved fcc-AI is subjected 
to pressures -60-100 GPa (see Methods for details). 



Indexing of the |iXRD pattern revealed that the observed Bragg 
peaks originate from a material having bec unit cells with a lattice 
constant a = 2.866 A and space group Im-3m (229). Such unit cells 
have only two unique atoms at atomic coordinates of (0, 0, 0) and 
(1/2,1/2,1/2). Hence, the two shortest interatomic distances are 
2.482 A (between the corner and centre atoms) and 2.866 A (between 
the corner atoms). As microexplosions were confined in the bulk 
of crystalline A1 2 0 3 , the only constituents that can contribute to 
the observed Bragg peaks are Al and O. The shorter interatomic 
distance agrees well with a double atomic radius of Al of 1.25 A 
(ref. 36), thus suggesting Al-Al bond. Because the atomic radius of 
oxygen is much smaller, 0.6 A, the length of Al-O bond would be 
about 1.85 A. For instance, bonds of similar length can be found 
in a crystalline A1 2 0 3 lattice, where the two shortest interatomic dis- 
tances are 1.852 A and 1.975 A (http://www.icdd.com/). In our bec 
lattice, the presence of Al-O bond would result in a jUXRD pattern 
drastically different from the one observed. Hence, oxygen cannot 
form bonds with aluminium in the identified bec structure with 
a = 2.866 A, and does not contribute to the experimental /iXRD 
pattern shown in Figure 2c. Consequently, both atoms in the bec 
unit cell are Al, and the observed /iXRD pattern is due to nano- 
crystallites of pure aluminium. The /iXRD provides a non-destruc- 
tive verification of the new bcc-AI phase which is stabilized in a 
compressed enclosure, as estimated by an observed co-existence of 
bec- as well as a -60 GPa-compressed fcc-AI phase as deduced from 
a shift of /iXRD peaks shown in Figure 4 (Methods). 

Spatial separation of ions. The observed formation of crystalline 
aluminium out of aluminium oxide is experimental evidence of a 
quite unusual and unexpected phenomenon, namely spatial sepa- 
ration of oxygen and aluminium ions in conditions of complete 
confinement and preserved stoichiometry. We propose the follow- 
ing explanation for this effect. Ultrashort laser pulses with intensity 
above the optical breakdown threshold break bonds of A1 2 0 3 and 
ionize the constituent atoms into plasma, which exists until the mate- 
rial cools down below the thermal ionization threshold. While this 
plasma exists, the Al and O ions diffuse and scatter at different rates, 
allowing spatial separation of the sapphire constituents (Fig. 5). 

The Rutherford cross-section, cr 12 °c Z^Z^fi ~ 2 (v x -v 2 )~ 4 gov- 
erns the Coulomb scattering of ions in plasma; here Z 1 and 
M M 

Z 2 are ion charges, /i = — - — — is the reduced mass of the 
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distance between Al- and O-ions can be expressed in the follow- 
ing form: 



Figure 5 | Synthesis of materials via spatial separation of elements 
in the plasma state. Schematic explanation of spatial and temporal 
transformations in the volume of ionized sapphire affected by a 
microexplosion, leading to spatial separation of Al and 0 ions and 
quenching of the bcc-AI. Horizontal and vertical axes represent time and 
distance from the centre of microexplosion, respectively. 



interacting particles, and (v 1 -v 2 ) is their relative velocity 37 . The 
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increases with the interaction energy. Hence, the mean free path 
becomes larger than the inter- atomic distance at a certain threshold 
energy level, around lOeV in our experiments. The difference bet- 
ween the scattering and the diffusion length of ions with differ- 
ent masses becomes apparent above the threshold. This effect is 
especially pronounced in collisions of particles with significantly 
different masses. The maximum difference occurs in the case of 
electron-ion collisions, where the smallest mass is equal to the mass 
of a single electron 37 . The difference in the diffusion length between 
oxygen and aluminium ions A = L° d - L A j leads to spatial separation 
that depends on the ion charge, mass and the rate of electron-to-ion 
energy transfer as estimated below. 

In a neutral solid, spatial separation of O and Al is physically 
impossible even when the molecular bonds are broken at elevated 
temperature, because the diffusion is insufficient to move atoms 
for a distance larger than the inter-atomic spacing. The major rea- 
son is that the cross-section for atomic collisions of neutral atoms 
is temperature independent. By contrast, collisions of ions in hot 
plasma are governed by the Rutherford cross-section, which is 
inversely proportional to the squared relative energy and reduced 
mass of the interacting particles 37 . The electron-to-ion energy 
transfer time in plasma is directly proportional to the ion's mass: 
C = M ion( m e v 7 m J> nere v 7 m is tne electron-phonon momentum 
exchange rate. Thus, the difference in mass becomes important at 
the early stage of plasma formation: the lighter O ions are heated 
more, gain more energy, and diffuse faster than the heavier Al 
ions. As a result, a unique non-equilibrium three-temperature 
state is established. The difference between the diffusion lengths 
A = L° d -L J d J leads to spatial separation of Al- and O-ions. The 

and diffusion velocity D d = / mpf v ion /3 



diffusion length L d = ^ 
is expressed through the mean free path / mfp and velocity v ion 
of the ion. Diffusion continues during the time t d until the plasma 
temperature drops below a certain level. Thus, the separation 
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(1) 



At a fixed temperature, the differences in the ions' velocity, mean 
free path, and diffusion length are related to the difference in their 
mass. The ratio of oxygen to aluminium ion velocity is expressed 
through the difference in mass as: v 0 /v A1 = (M A1 /M 0 ) 1/2 = 1.3, con- 
sequently, the spatial separation can be expressed as follows: 
A ~ 0.33L°; we assume here that each ion diffuses in plasma with 
the average ion mass M av = 20 a.u. and jl Al0 = M AW MJ(M AW + M av ). 
One can see that the diffusion length of oxygen should exceed 
100 nm to form the tens-of-nm sized crystallites observed in the 
experiments. 

Discussion 

Let us now estimate the diffusion coefficient of oxygen and the 
diffusion time t d , which define the diffusion length. The Coulomb 
collision cross-section for the temperature range of 10-50 eV only 
slightly exceeds the atomic cross-section. It is around <7 R ~10~ 15 cm 2 . 
Therefore, the mean free path in solid density plasma is comparable 
to the interatomic distance. The lifetime of the plasma limits the ion 
separation period. Thermal ionization by electron impact continues 
the ionization process after the pulse almost at the same level as dur- 
ing the pulse, because the plasma temperature is much higher than 
the ionisation potential 29 . Therefore, plasma exists until cooling 
processes due to energy dissipation in the expanding shock wave 
and heat conduction reduce the temperature below the threshold 
of-lOeV. 

The process of void formation by rarefaction wave maintains the 
high temperature of ions in the central parts of the laser- affected 
volume. The density decrease in the central part of a spherical 
explosion is accompanied by an increase in temperature, keeping 
the pressure almost evenly distributed over the space of the central 
region 37 . The ion's mean free path is inversely proportional to the 
plasma density. According to hydrodynamic simulations of micro - 
explosions 29 , mass density at the centre is around 0.1 gem -3 at 0.9 ns 
after the microexplosion. Hence, the mean free path of oxygen ions 
moving away from the rarefied central area is within several nanom- 
eters, and their diffusion rate is around D°~0.1 cm 2 s _1 . Taking the 
cooling time t d ~l ns in accordance with computer simulations 29 , the 
oxygen diffusion length is Lj~100nm. Thus, separation of oxygen 
and aluminum ions in hot and dense plasma due to the difference in 
the diffusion rates appears to be a viable process (see equation (1)) 
responsible for the formation of superdense bcc-AI crystals. 

To conclude, this study shows that a femtosecond laser- induced 
microexplosion confined inside sapphire leads towards synthesis of 
a new and previously unobserved dense bcc-AI phase, which sur- 
vives in a compressed state after fast quenching. The route of synthe- 
sis is via spatial separation of Al and O ions in short-lived hot non- 
equilibrium plasma of solid state density. We have demonstrated 
that high energy density produced in a simple tabletop experiment 
makes it possible to form an exotic high- density material phase, 
which could not be produced by other means. The mechanism 
responsible for the formation of the new phase of material can occur 
in a wide range of other materials that could be similarly subjected 
to extreme pressure and temperature conditions and converted into 
spatially confined plasma. As ultrafast lasers can operate at very 
high repetition rates within the MHz- range and deliver a high pulse 
energy 38 , laser-induced micro explosions offer a new strategy for the 
synthesis of high-pressure phases of materials in useful amounts in 
table-top laboratory conditions. 
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Methods 

Laser modification of sapphire. The laser source was a Hurricane X (Spectra 
Physics) regeneratively amplified Ti:Sapphire femtosecond laser system, operating 
at a wavelength of 800 nm and producing pulses with a temporal length of about 
150 fs. The pulses were focused in the sample using an oil-immersion objective 
(Olympus) with numerical aperture NA= 1.4. The pulse energy in the focal spot 
was 135±6nJ, corresponding to a peak intensity of I 0 = 487 TW cm -2 . The sample, 
an 80 um-thick sapphire slab (Shinkosha), was mounted on a high-precision 
three-dimensional translation stage controlled by piezoelectric transducers (Physik 
Instrumente). The stage was translated between consecutive laser shots along a 
predefined trajectory, such that an ordered pattern of void-structures was produced 
in the slab. The focusing depth below the surface of the sapphire sample was within 
the range 20-60 um. Precompensation of the pulse group-velocity dispersion in 
the microscope lens and other optical elements along the optical path was done to 
ensure the shortest pulse duration at the focus. 

X-ray diffraction experiments and data processing. /iXRD characterization 
of the laser-processed region was carried out by angle-dispersive X-ray microdif- 
fraction at HPCAT beamline 16-BM-D of the Advanced Photon Source, Argonne 
National Laboratory. The monochromatic X-ray beam with energy of 33keV 
(A = 0.37571 A) delivered by a Si (1 1 1) monochromator was focused by a pair 
of Kirkpatrick-Baez mirrors down to an area of (5.0x10.0) \im 2 (widthxheight), 
estimated at FWHM. X-ray diffraction patterns were collected in the transmis- 
sion mode using a Mar345 Image Plate with typical exposure time per image of 
50-100 s. The sample-to-detector distance and detector tilts were calibrated with a 
Ce0 2 (NIST) standard. All two-dimensional /iXRD images were integrated to one- 
dimensional diffraction patterns using FIT2D software. Theoretical powder X-ray 
diffraction patterns expected for fcc-Al, bcc-Al, sapphire and high-pressure A10 x 
phases were calculated using the Materials Studio package (Accelrys). 

XRD pattern analysis, structure determination and Rietveld refinement was 
performed using X'Pert HighScore Plus (PANalytical) and Materials Studio 
(Accelrys) software. By indexing the X-ray diffraction pattern, we unconditionally 
identified the size and symmetry of the unit cell. The ability to observe seventh 
peak in the jiiXRD pattern enabled us to avoid ambiguity and distinguish simple 
cubic and bcc structures. The unit cell obtained is body-centred-cubic with the 
unit cell parameter a = 2.866 A and space group Im-3m(229). The peak widths and 
intensities were fitted using Rietveld refinement. Peak shapes were fitted using 
crystallite size and strain as fitting parameters. This resulted in a crystallite size of 
17.9 nm and strain of 0.26%. The resulting Rietveld refinement profile of the bcc-Al 
phase is shown in Figure 3. The last Bragg peak intensity was not included in the 
fit, because only a fraction of the diffraction ring was acquired in the /iXRD image. 
The R- indices of the criteria of fit are R^ p = 0.08 and R^ = 0.05, which represent a 
high-quality fit for compounds with high symmetry. 

The major jtiXRD peaks ride on a weak broad background, which also features 
many small peaks. This background is emphasized in the plot in Figure 4. The 
broad background feature at small diffraction angles (Fig. 4) originates from 
amorphous sapphire created by microexplosions, whereas the sharp not-Al-related 
peaks most likely originate from oxygen-rich A10 x phases as well as damaged and 
misoriented A1 2 0 3 regions. The fcc-Al phase was also present in the laser-affected 
region. jiiXRD spectra revealed the presence of compressed fcc-Al diffraction peaks 
in several void-and-shell structures (see Fig. 3). The set of the diffraction peaks 
corresponding to fcc-Al and bcc-Al are compared in Figure 3b. As it can be seen, 
the fcc-Al pattern belongs to the space group Fm-3m (225) (http://www.icdd. 
com/), the same as that for fcc-Al at the ambient conditions. However, the atomic 
volume of observed fcc-Al, V= 1 1.51 A 3 , is 1.44 times smaller than the volume of 
fcc-Al in the ambient conditions, V 0 = 16.62 A 3 . Thus, in our experiments fcc-Al is 
compressed practically to the same density as the newly born bcc-Al. The Young 
modulus of aluminium at room temperature comprises 70 GPa. An increase in 
density would result in a rise of the Young modulus up to 100 GPa. We estimate 
cautiously that residual pressure that will keep compressed fcc-Al inside a sapphire 
matrix is in the range 60-100 GPa (ref. 22). Therefore, one may expect residual 
strain to exist in the laser-compressed shell where fcc-Al and high-density bcc-Al 
were observed. 

Hume-Rothery analysis of Al and O solid solution. The XRD profile compa- 
rable to that of bcc-Al can also be generated by a bcc unit cell of a substitutional 
alloy of Al and O with the same lattice constant a = 2.866 A. In the case of Al and 
O occupying distinct bcc lattice sites, owing to dissimilar scattering factors of 
Al and O, XRD pattern will exhibit Bragg peaks such as (111), (201), (122), (003), 
(113) and (203) within 26= 10-30° that are otherwise forbidden in bcc crystal 
structure. As these peaks are not present in the experimental XRD pattern shown 
in Figure 2, we therefore exclude the existence of ordered Al-O alloy phase in our 
sample. The Bragg peak positions and relative intensities will match those of pure 
bcc-Al only, if both atomic sites of the bcc unit cell are randomly occupied by 
Al and O and with the same Al/ O ratio, similarly to the case of Ce and Al 39 . The 
formation of such solid solution, however, has to satisfy empirical Hume-Rothery 
rules, which state that the atomic radii of the solute and the solvent must not 
differ by more than 15% and that their electronegativity should be within 0.4 of 
each other 34 . 



To calculate atomic radii of Al and O we use the phenomenological equation of 
state (EoS) by Vinet et al. 40 : 

P(x) = 3B 0 [((1 - *)) I x 2 ] exp [1.5(fl5 " D(l " *)L ( 2 ) 

where B 0 is the bulk modulus, B\ is pressure derivative of the bulk modulus, 
x = 3jV I V 0 » is tne initial and Vis the high-pressure phase volumes, respec- 
tively. We calculate atomic radii of Al and O at two pressures: 100 GPa — the 
estimated maximum remnant pressure of bcc-Al in A1 2 0 3 matrix after the laser 
pulse and 275 GPa— the theoretically predicted pressure for bcc-Al formation 41 . 
Following the EoS obtained by Akahama et al. 22 , we estimate atomic radii of 
bcc-Al as 1.189 and 1.087 A, at pressures 100 and 275 GPa, respectively. To form 
solid solution with Al, oxygen has to be in a monatomic phase with atomic radii 
within 15% of that of Al. However, up to this date, no monatomic phase of solid 
oxygen has been experimentally observed up to pressures of 1 16 GPa (refs 42,43) 
and therefore experimental EoS for monatomic oxygen does not exist. Moreover, 
it has been predicted by Oda et al. 44 that molecular (diatomic) oxygen phase is 
stable over a wide range of pressures 100-2,000 GPa. The high stability of diatomic 
oxygen up to very high pressures makes it very unlikely that substitutional Al-O 
alloy can be formed. 

Nevertheless, some theoretical work has been done to predict monatomic solid 
oxygen phase and estimate its atomic radius at high pressures 44 " 46 . In the work by 
Anderson et al. 45 EoS of Fe oxides, corrected to remove the effects of Fe, was used. 
They estimated an atomic radii of monatomic oxygen at pressures 136, 330 and 
360 GPa to be 0.94, 0.89 and 0.88 A, respectively 45 . Because elements in many groups 
of the periodic table under high pressure exhibit similar sequence of structural 
phase transitions, further theoretical predictions utilized the fact that other group 
Via elements, S, Se and Te, do form monatomic phases at high pressures 46-51 such 
as /3-Po rhombohedral phase, which can be considered as simple cubic (sc) phase 
that is deformed along [111] direction. According to the band structure calculations 
performed by Suzuki et al. 46 , the most stable /3-Po unit cell shape of monatomic 
solid oxygen is when cla= 1.04, where c and a lattice constants refer to the hexago- 
nal setting 46 . Using this result and applying EoS obtained for molecular oxygen 42 , we 
estimate radius of monatomic solid oxygen in /3-Po crystal structure as 0.906 and 

0. 827 A for pressures 100 and 275 GPa, respectively. These values are in good agree- 
ment with the atomic radii reported by Anderson et al. 45 The atomic radii difference 
between Al (solvent) and O (solute) at pressures 100 and 275 GPa are 24% and do 
not satisfy HR rules for solid solution formation. Following the trend of group Via 
elements, one can predict a formation of bcc- O crystal structure at even higher 
pressures 22 ' 50 . Monatomic oxygen in such structure will display atomic radii of 0.986 
and 0.902 A at pressures 100 and 275 GPa, respectively. For the bcc-Al and bcc-O 
structures, the difference in atomic radii is 17%, too great to satisfy HR rules. 

The electronegativity can be calculated as EN= 0.31x(Z+ l)/r c + 0.5 (ref. 52), 
where Z is the valence electron numbers (3 and 6 for Al and O, respectively) and 
r c [A] is the covalent radius. Assuming bcc-Al as a covalently bonded material, 
the estimated electronegativities for Al are 1.543 and 1.641 at pressures 100 and 
275 GPa, respectively. Considering that monatomic oxygen in /3-Po and bcc crystal 
structures is also covalently bonded, we can get the following electronegativites 
for O: /3-Po structure— 2.896 and 3.125, bcc structure— 2.700 and 2.906 at 100 and 
275 GPa pressures, respectively. The differences in electronegativities between Al 
and O are much larger than 0.4 for all estimated structures and therefore do not 
satisfy HR rules for solid solution formation either. 

On the basis of XRD data analysis and estimates of atomic radii and electro- 
negativities, we conclude that solid solution of Al and O cannot be formed. 
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